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1. Introduction 

Quantum decoherence is at the heart of both the foundations and applications of 
quantum physics. Cavity quantum electrodynamics (QED) systems, operating in the 
strong coupling regime, have proven to be excellent for the studies of the entangled atom- 
photon state [1]. The experiment involving the Rydberg atoms in a high Q microwave 
cavity have opened the way to the studies of the decoherence dynamics in a mesoscopic 
system [2] . Theoretically, the simplest model that captures the physics of such a hybrid 
system is the Jaynes-Cummings (JC) model [3|,|1]. It is one of the few exactly solvable 
models in quantum optics and predicts several interesting effects such as the vacuum 
field Rabi oscillations O El [7] , collapses and revivals of Rabi oscillations in the coherent 
field |8]. Remarkably, it was noticed that the atom is to a good approximation in a pure 
state in the middle between the collapse and revival [HI [10] . 

On the other hand, being related to the quantum measurement theory and 
the quantum decoherence problems, the influence induced by the spin bath on the 
decoherence dynamics of a central system have also attracted much attention [llj. It 
was shown that the decoherence induced by coupling a system with an environment may 
display universal features: the decay of quantum coherences in the system is Gaussian 
for the specific initial environment state [121 Il3l HI] 

In this paper, we revisit the decoherence of the atomic state in the resonant JC 
model with the field initially being in a coherent state and elaborate the dynamic 
evolution of the purity. Closed analytical expressions for the purity of the central 2- 
level atoms are obtained. We observe that the similar behavior as that of the spin bath 
occurs in such a hybrid system. Both the first collapse and the amplitude of subsequent 
oscillation exhibit a Gaussian decay behavior. 

The paper is organized as follows. In Sec. II, we introduce the main properties 
of JC model. In Sec. Ill, we evaluate the purity of the central 2-level atom. Section 
IV presents our summary and conclusion. In the Appendix we derive the main formula 
needed to evaluate various sums used in the text. 

2. The JC model 

Starting point of the analysis is the JC model, which consists of a single atom coupled to 
a single mode cavity. The two possible states of the atom are the ground state \g), and 
its excited state |e). The model and the subject we discussed are the same as the one 
previously studied by Gea-Banacloche [9l HO], who showed that the atom is to a good 
approximation in a pure state in the middle between the collapse time tc and revival 
time tr- In the following, we will not only reproduce the results of [9l [10] but also show 
that this remarkable phenomenon is a part of the dynamical process for a special case. 
The model Hamiltonian at the resonance has the form 




(1) 
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= = |e) {gl , (^z = |e) (e| - |y) (^| 

where is the creation operators of photon with frequency a;, a;a is the atomic transition 
frequency, and A is the cavity-atom couphng constant. The aim here is to study the 
dynamics of a given initial state. A general initial state of the system has the form 

oo 

\^{Q)) = {C,\g) + Ce\e))^Cn\n) (2) 

n=0 

where |Cep + |Cg|^ = 1. Of central importance is the excitation number 

^^=^(^z + a^a + i (3) 

is a conserved quantity, i.e., [J\f,H] — 0, which makes it easy to diagonalize the 
Hamiltonian, since the atom-field eigenspaces are only two-dimensional. It also makes 
the dynamics of states involving several subspaccs simple. Nevertheless, the dynamics 
of states that have many significant energy-state components can show considerable 
complexity. 

Introducing a unitary transformation 

7e(^,0) =e^K/^+'^"^") (4) 

which generate the phases 9 and on the atomic and cavity states 

n {9, ct>) (C, \g) + Ce |e)) \n) = e^^ (e-^/^C, |^) + ^'I^C, |e)) \n) , (5) 
we have TZHTV = H, where 



H — X (cr^a + (7_a^j -|- -uJa<^z + ojo^a (6) 

— A {dj^a + a"_a^^ + -uJaC^z + ooa^a 

with a = e~^^'^~^^a, 5"+ = e~'*^'^~^V+ and = a^. It indicates that Hamiltonians H and 
H share the same eigenfunctions by transformation of the basis {\n)} — )■ |e'"''^|n)| 
or \g) — )■ e~'^/^ \g) and |e) — )■ e'^^^ |e). Remarkably, in the case of ^ = 0, we 
have TZHTZ^ = H = H, which shows the invariancc of the Hamiltonian under the 
transformation TZ [9, 9). We will show that, when dealing with the coherent cavity state, 
this feature leads to an interesting and important phenomenon. We will demonstrate 
the strong dependence of the dynamics of the atomic purity on the relative phase of the 
atom and the cavity field. 

We shall only consider the resonant case of = Ua- Then at time t, state |^ (0)) 
evolves to 



oo 

I* (0) = E [CnCg cos (V^At) - iCn-iCe siu (v^Ai)] \g, n) (7) 

n=l 

oo 

+ E [^nCe COS (VrT+TXt^ — \Cn+iCg siu (^V^T- + lAt^ |e, n)]. 



n=l 
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We concern the reduced density matrix of the atom, which has the form 

a b 
b* l-a 



P.W-f" M (8) 



where 

oo 

J2 [\CnC/ cos' (v^At) + \a,-iCf sin' (^^Xt) (9) 



n=l 



-2Im {CnC^_^CgC*) COS (v^At) sin (v^At 

oo 

[\Cnf CgCl COS (v^At) COS (y n + lAt I 

n=\ 

+Cn-iC^+iCeC* sin (VnAt) sin (V^ + IXt 



+iCnC*_^_i \ Cgf COS (^v^At^ sin (^y/n + lAt^ 

— iC„_iC* |Ce|^ sin (^v^-^^) cos (^Vn + lAt^ 
As a measure of the degree of coherence, the purity the atom can be expressed as 

Pit) = Tt (pI) = a' + (1 - a)' + 2\b\\ (10) 
where Tr(...) denotes the trace on the cavity field. 

3. Decoherence of a two-level atom 

With the time evolution of the reduced density matrix of the atom, we can investigate 
the dynamical behavior of the atom, which has been employed to calculate the inversion 
and the purity for the case of initial coherent state. The initial state has the form 



\^m = iCg\g) + a\e))\a) (11) 



\a) = e-l"l /' ^ ^ |n) , a = e-'^^Vfi. (12) 
n=o vn! 



where 



It is has been found that the initial Rabi-oscillations concerning the probability of being 
in a given atomic state decay on a timescale called the collapse time, tc = 2/X, but then 
revive after a much longer time, tr = 27Ty/E/X [151 US]- 

Here we discuss the time dependence of the atomic purity in a long time scale. For 
\af (or n) ^ 1, we have 

Cn = e-l"l'/'^ = ^a-i ^ e-^^C„_i. (13) 



Remarks on flat condition. In the following we take = for the sake of simplicity, 
since factor e""*^ can be mapped on the atomic state by {Cg, Ce) — )■ {e'^'^^'^Cg + e"^/'Ce^ 
according to our previous analysis. Then in the following derivation, we simply consider 
the coefficients Cg and Ce as complex numbers. On the other hand, for sufficiently large 
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value of n, the Poisson distribution is an approximation to the normal (or Gaussian) 
distribution, i.e., 



1 



exp (—n) — ^ 
nl 



V2 



= exp 



nn 



[n 



-\2n 
n) 



2n 



(14) 



On the other hand, the Poissonian function peaks sharply around n. Then for a 
nontrivial function F(y/n, 1/y/n), one can take the approximation 



V2 



irn 



: exp 



[n 



n) 



2n 



n 



n 



n 



2 2Vfi'2y/fi 



=f)-(15) 



Furthermore, in the limit of n ^ 1 the summation over n can be done exactly by virtue 
of Euler-Maclaurin formula. In the Appendix we derive the main formula needed to 
evaluate various sums used in the text. Accordingly, we obtain the analytic form for 
reduced density matrix elements are 
1 



+ 



cos 

V tc J 



1^ ^ I • n ■ ( 4:y/nt\ 
+ \CgCe\ smt^sm I — — — j 



/ 2t'' 



|CgCe| COS COS 



I + i- sin 



exp 



'Sn tl 



+ i 



^ (iC'gl - \Ce 



sm 



+ \CgCe \ sin 6 cos 



2£^' 

Tl 



From equation ffTOj) . we have 



Pit) 



1 

2 ' 
+ 



|CgCe|^ COS^ COS^ 



\Cg\'-\CX) +4|C,Ce| sin^^ 



TTt 



exp 



T^t^ 

An tl^ 



exp — 



f2 



'Tlt\ 

+ sin I — 



Cg\ — \Ce 



+2 \Cg\ \Ce \ sin 9 cos 



. fA^/f^t\ 

sm 

V tc J 



mti 



exp I - 



2f' 

Tl 



(16) 



(17) 



We can see that, at small time region t <^ tr, the term containing exp (— 4t^/t^) is a 
transient state and is dominant at the beginning of the evolution. This Gaussian decay 
behavior is independent of the mean photon number n. After the transient relaxation, 
the purity of the atom exhibits oscillating Gaussian dependence on the time. The width 
of Gaussian function is independent of the values of Cg and Ce- The initial atomic state 
determines the amplitude of the oscillation. Both the first collapse and the amplitude of 
subsequent oscillation exhibit a Gaussian decay behavior. Using the analytic expressions 
for the purity amplitude equation ( IT71) . we estimate the period of the oscillation to be 
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Figure 1. (Color online) Plots of the purity obtained by the exact diagonalization 
(solid line) and the approximate analytical expression equation (|17|) (empty circle) at 
various relative phases 6, for the initial atomic state with \Cg\ — \C'e\ and the initial 
coherent state of photons with average number n = 400 and n = 16. 



the same as that of the Rabi oscillation. Obviously, purity dynamics depends on the 
parameter of the system as well as the initial state. In this work we show that the 
relative phase between the initial atom and cavity field has far more important influence 
on the purity dynamics. Let us consider two interesting special cases: |C„Ce| = and 



\Ce\- In first case, = ( or Ce = 0), we have 



Prr 



.W = o + oSm — exp - — - 



4n tl 



1 / 4t2^ 



'in — — 



-K^t^ 



exp 



2e 
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Figure 2. (Color online) The same as figure [T] but in the small time scale and the 
case oi 9 — 7t/2. 



i.e., the amplitude of the oscillations becomes maximum. In second case, Cg 
we have 

1 1 



Peit) = 7; + 7; 



1 — sin^ 9 cos^ 



+ -sm ^exp I 
+ sm j sm tl cos 




(19) 




It shows that after the transient process, the amplitude of the oscillations only depends 
on the relative phase 6. We also note that Pe{t) behaves as the purity for various values 
of the coefficients Cg = e"'^/^ \ Cg\ and Cg = e^^/^ |Ce| by simply replacing sin^ 6' in the 
equation fl20l) with 1 — 4 \CgCe cos . Thus in the following numerical simulations, we 
only demonstrate the case of \Cg\ = \Ce\ for simplicity. 
We note that in the case of ^ = 



,11 / nH^' 



(20) 



which corresponds to the envelope of the pattern P{t) for arbitrary initial atomic state. 

The above analysis shows two important characteristics of the decoherence 
dynamics. At first, the decoherence occurs dramatically at the very beginning for an 
arbitrary initial atomic state except the case of Cg = Ce- Secondly, after the transient 
decoherence, the amplitude of the oscillating purity strongly depends on the initial 
phase difference between the atom and the field. When the initial atom and the field 
are in-phase or opposite-phase {(p — 6 = 0, tt), the atom has a relatively long coherent 
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Figure 3. (Color online) Plots of the purity for the top-hat initial field states with 
average number n — 400, D = 10, 20, 30, and 80, and the approximate analytical 
expression (empty circle) for the initial coherent state of the field with the same 
average number. It indicates that the atomic coherence decay differs from the Gaussian 
function except the special case. 



time. When they are orthogonal-phase {(f) — 6 = 7r/2, 37r/2), the atom acquires maximal 
oscillating amplitude of decoherence. 

In order to verify the above analysis some numerical simulations are performed. In 
figure [H we plot the equation ffTOj) for n = 400 and 16 cases with different values of 6 
and \Cg\ = \Ce\- As comparison, we also plot the equation f lT7|) accordingly. We can see 
that the analytical results match well with the simulation results, especially in large n 
case and during the first several periods of oscillation, figure [2] is the same as the plot 
in figure [U but for small time scale and 9 = 7r/2 to demonstrate the transient process 
explicitly. 

It is also worthwhile to mention that the Gaussian decay of the decoherence is the 
direct result of the coherent state environment. We note that the expression equation 
( !T7|) is obtained under the two conditions: (i) distribution function C„ is flat as equation 
( IT3ll : (ii) we can use the approximation equation ( IT5l) near the mean photon number 
n. The result for 6 = may promise important potential applications in quantum- 
information processing since Gaussian time dependence of the decoherence factor would 
suggest a different more frequent error correction than the exponential dependence. 
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4. Summary 

In conclusion, considering a system consisting of a two-level atom, initially prepared in 
a coherent superposition of two levels, interacting with a coherent state of the field, we 
show that the dynamics of the atomic purity are sensitive to the relative phase between 
the atom and the cavity field. We also observe that the purity of the atom exhibits 
oscillating Gaussian dependence on the time with a width independent of the initial 
atomic state. Our results may have a great potential for future applications in quantum 
optical device. 



Acknowledgment 



We acknowledge the support of the CNSF (Grant Nos. 10874091 and 2006CB921205). 

This appendix contains the formulas needed to evaluate various sums used in the 
paper. 

Sum 1. Consider the sum 

[n — nf , Xt 



Si 



J^exp 



According to equation (IT^ . we have 

oo 



2n 



+ 1 



2^^ 



V2 



Tcn 



n=l 



{n — n) I . At / 



2n 



+ 1 



From Euler-Maclaurin formula, we replace the sum by the integral 



Si 



1 



y/27mJinfty 

( XH^ 

exp 



exp 



(x — ny 



2n 



. Xt 
4a/^ 



(-1) 



dx 



. Xt \ 
32^2 + ^27i / 
where the Gaussian integral formula 



a\j2'n J -<x> 



exp 



2cr2 



— \ 2 



-i/3t 



d/3 = exp 



-i/3t 



(-1) 



(-2) 



(.3) 



(•4) 



has been used. 

Sum 2, 3. Taking the similar procedure one can calculate the following two sums 

\2- 



and 



oo 

V 27rn 



oo 

53 = E e^P 



V^n=l 



{n — ny 
2^ 



[n — ny 
2n 



cos \/n 



At) 



cos {^^/nXt^ sin {y/n + lAt) . 



(.5) 



(.6) 
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Actually, using the Gaussian integral formulae 



and 



we have 



and 



1 /■°< 



exp 



— \ 2 



ay/2TT J-oo 



So 



exp 



2(72 



2cr2 



cos {/3t) d/3 — cos (^/3t^ exp 



sin d/3 = sin ((3t) exp 



2+2' 



2+2' 



S. 



V2 



f 

J — I 



exp 



{x — nY 



1 ^ 1 

— oo 
2J.2" 



2n 



cos 



dx 



i + exp 



exp 



cos 



(x — nY 
2n 

{ynXt) 



COS 



n + ^ At 
VnJ 



dx 



2V2 



+ sin 



2 exp 



T^nJ-oo 

X 



exp 



[X — n) 



\2l 



2n 



xXt 



sm 



3Xt \ 



X 



,Xt+[^/f^+ —j= Xt 



X 



+ - exp 



nt'^ ( X 



2 \\/n A\fr^ ^ 

nt^ ( X \ 



sin 2\/nXt + 



jdx 

At ^ 
2^/n) 



2 l4v^, 



sm 



At ^ 



2+2' 



A^t 



2 exp 



sin(2V^At) + iexp(-Al_jsin(^^ 



■7) 



(.9) 



.10) 
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